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Abstract 


We show that it is possible to model two-body resonant interactions at low energy with a class of finite- 
range potentials based on the methods of lost and Kohn. These potentials are expressed in terms of the 
effective range tq and the s-wave scattering length a<j. We derive continuum solutions of these potentials. 
By writing V±{r) = Vo(^) + where the sign -i-(-) refers to positive(negative) scattering length, 

Vo{r) is of the form of Pdschl-Teller potential and is expressed as a power series of the small parameter 
e = (-^1 — 2ro/as)~^ — 1 when is large, we derive Green function of Vo(r). Using the Green function, 
solutions of V±{r) for |as| » tq can be obtained numerically by treating V^{r) as a perturbation. We de¬ 
scribe the threshold behavior of scattering phase shift for Vo(’’)- This study may be important for developing 
a better understanding of physics of strongly interacting ultracold atomic gases with tunable interactions. 
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I. INTRODUCTION 


Physics of interacting many-particle systems is deseribed, to a first approximation, in terms of a 
mean-field potential. At low energy, there is a well-known mean-field potential of zero range. This 
is known as Fermi’s eontaet potential or pseudopotential yj], usually expressed as a delta funetion 
with the potential strength being proportional to s-wave seattering length a* only. For sueh a 
delta-type potential, s-wave seattering wave funetion at zero energy beeomes singular at r = 0, 
where r is the separation between two partieles. To eireumvent this problem, eontaet interaetion is 
replaeed with a suitable bound^ eondition at r = 0 in the Sehrodinger equation resulting in the 
regularized pseudopotential H-lSD 

Airfi^n., r) 

( 1 ) 


d 

1/ r = -d r —r 

^ ^ 2/i ^ ’ dr 


where /r is the redueed mass of the two partieles. The zero-range potential method is applieable 
to low temperature dilute atomie or moleeular gases whose two-body interaetions obey Wigner’s 
threshold laws h. For two spherieally symmetrie atoms, Wigner’s threshold laws dietate that 
beeomes an energy-independent parameter. The s-wave seattering amplitude then reduees to the 
form /o(A;) ~ whieh is /o(A;) ~ for k\as\ « 1. However, when a quasibound or 

bound state exists near the threshold of a potential or a seattering resonanee oeeurs at low energy, 
as beeomes large leading to the breakdown of the eondition k\ as\ « 1. Wigner’s threshold laws 
do not apply when a resonanee oeeurs near zero energy. As a eonsequenee, Fermi’s pseudopoten¬ 
tial approaeh beeomes inadequate to deseribe eollision or many-body physies near resonanees. 

The main purpose of this paper is to seek a finite-range model potential that ean deseribe 
resonant interaetions. Resonanees and physieal phenomena related to resonant interaetions are 
ubiquitous in physies and ehemistry. In ultraeold atoms, Feshbaeh resonanees ean be indueed 
by external magnetie [lZL ISD or optieal 11914 1111 or both magnetie and optieal fields [Il2n . In eur- 
rent literature, many-body phenomena in ultraeold atomie gases near a Feshbaeh resonanee are 
by and large deseribed in terms of eontaet potential [Il3n . However, the effeetive range and en¬ 
ergy dependenee of many-body phenomena ean hardly be ignored when a resonanee oeeurs. It 
is theoretieally shown that for » i?vdw when is mueh larger than the eharaeteristie length 
seale i?vdw = 12 of van der Waals potential, the effeetive range tq is about three 


times i?vdw 1114, 


1511 . The values of R^dw for different alkali atoms used in eold atom experiments 


are tabulated in Ref. [SD. For » i?vdw> the effeetive range tq of a potential with van der 


Waals dispersion at long separation will be several tens of Bohr radius (oq). This means tq is large 
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compared to typical molecular equilibrium positions. It has been recently shown that the range 
and energy dependence of magnetic Feshbach resonances of ultracold atoms are quite important 


m 


ITII . Particularly in case of narrow Feshbach resonances, tq can be very large, even of the or¬ 


der of hundred or thousand Oq. Currently attempts are being made to construct an improved model 
potential by resorting to a toy model lIlSll or contact potentials with energy-dependent scattering 
length [fl^ . The fact that a contact interaction can not accurately describe resonant interactions 
calls for the formulation of a better model potential that can take into account both the finite range 
and the energy dependence of scattering amplitudes. 


Here we show that one can deduce, based on the methods of lost and Kohn OOl 


2111 . a class 


of finite-range model potentials that can account accurately for s-wave resonances at low energy. 
We demonstrate that, for large a^, the model potentials reduce to Poschl-Teller form ll22ll that can 
admit analytical solutions. Using these analytical solutions, we construct Green function that can 
be used for solving the full potential numerically. The primary physical motivations behind this 
work is to establish the connections of strongly or unitrity-limted interactions with the exactly 
solvable Poschl-Teller potentials. This study showing that the unitarity-limited interactions can be 
treated with exact analytical solutions may be important for gaining new insight into the physics 
of strongly interacting ultracold atomic gases. There is another motivation behind undertaking this 
work. In a seminal paper Butsch et al. [l23n have presented exact solutions of two identical cold 
atoms interacting via regularized contact potential in one, two and three dimensional harmonic 
traps. However, as expected, these exact solutions will not apply to strongly interacting atoms 
in traps due to finite-range effects. It is therefore worth seeking exact solutions of two atoms 
interacting via the finite-range model potentials with large scattering length in a harmonic trap. 
Such solutions will lead to new insight into physics of strongly interacting systems in confined 
or low-dimensional space. In this paper we do not make any attempt to find solutions of two 
trapped atoms interacting with finite-range potentials. However, the finite-range potentials and 
their exact solutions presented in this paper may serve as a precursor towards generalizing the 
results of Butsch et al. 11231] for finite-range potentials. 

This paper is organized in the following way. In the next section, we present our proposed 
finite-range potentials for resonant interactions and discuss how they are obtained. In Sec. 3, we 
present continuum and bound state solutions of the potentials. Analytical and numerical results 
are discussed in Sec. 4. The paper is concluded in Sec. 5. 
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II. FINITE-RANGE MODEL POTENTIALS 


The model potentials we eonsider are some variants of Bargmann potentials [|24ll and derived 


using effeetive range expansion of phase shift. The eonneetions between phase shifts and poten¬ 
tials were first rigorously established by Bargmann OSll . The early works on the method deriva¬ 
tions of finite-range potential from the phase shift data were earried out by GeTfand and Levitan 


ll26ll . lost and Kohn Od. 


2111 . and many others. Here we follow the method of lost and Kohn. A 


partieularly simple model potential results in when the phase shift 5o{k) is given by effeetive range 
expansion 


cot So{k) = 


ka 


1 1 , 

+ T^^ok + 


( 2 ) 


at low energy under the eonditions kro < 1 and tq « |as|. Jost and Kohn ll2flll showed that, for 
negative a^, the model potential takes the form 

Ah? a(3'^exp{-2(3r/ro) 


Vjr) = 


(3) 


fiVQ [a + exp(-2/9r/ro)]2 
where a = a/1 - 2ro/as, /3 = 1 -f a and /i is the redueed mass. This is valid for /o(-E')| < 7r/2 
in the limit E ^ Q. When a seattering resonanee oeeurs, k\ a,| » 1 sueh that if one negleets the 
effeets of effeetive range then 5o ~ 7r/2. This means that the S-matrix element exp[2i(5o] — — 1 
leading to unitarity-limited interaetions. We define unitarity regime by —1 « {kas)~^ « 1. 

When as is positive and large, the potential ean support one near-zero energy bound state. 
Therefore, to obtain a model potential for positive a* from effeetive range expansion, one needs to 
ineorporate the binding energy of the bound state. Then the potential beeomes a three-parameter 
potential. This potential is not unique unless an additional parameter eorresponding to the bound 
state is taken into aeeount. As shown by Jost and Kohn 112ill , the s-wave binding energy Eb = 
—h?KQ/2fi (where kq > 0) ean be parametrized, under effeetive range expansion, by introdueing 
a parameter A to express kq in the following form 


Kq = 


To L 


1 + Vl - 2ro/c 


1 + A 
J 1 - A 


(4) 


It is bounded by —1 < A < 1 for > 2ro and tq > 0. Using these three parameters a^, tq and A, 
an expression for the model potential is given in Eq. (2.29) of Ref. [l2lll . where a is bounded by 
0 < a < 1. Now, if we make a ehoiee A = —\/l — 2rQ/as the Eq. (2.29) of Ref. Elll reduees to 
a two-parameter potential of the form 

4/12 a/32exp(-2/3r/ro) 


U+(r) = 


/ir^ [1 -f aexp(-2/3r/ro)]^ 


(5) 
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This choice of A corresponds to the binding energy Eun — (2/ia^) for ‘IrQjas « 1. 

It is easy to notiee that in the limit —)■ ±cx), both the potential of Eqs. dS]) dS]) reduees to the 

form 

4^2 


Ko = 


/irQ cosh^(2r/ro) 


( 6 ) 


This form of the potential has been employed by Carlson et al. [l27h for quantum Monte Carlo 


simulation of a homogeneous superfluid Fermi gas with infinte negative seattering length. Shea et 
al. [12811 have ealeulated the energy speetrum of two harmonieally trapped atoms interaeting via the 
potential of Eq. d6l), showing that bound state speetrum remains almost the same as for zero-range 
pseudopotential ll23ll if tq is mueh less than the eharaeteristie length seale Iho = \/h/fiojho of the 
harmonie oseillator trap with frequeney ujho- 

The potentials of Eqs. d3]) and dS]) ean be written in the forms V±(r) = Vo{r) + V^\r) where 

Vo{r) = -^-(7) 


yW COsh"^ [ kt ] 


and 


rj"’ = ^"o5^(-ir(n + i) 


n=l 


1 + exp(±2fi;r) 


( 8 ) 


V^ir) is in the form of Posehl-Teller potential of seeond kind. Posehl-Teller potentials and their 
different variants are well-known in quantum meehanies 1129113011 as exaetly solvable potentials in 
one dimension. In three dimension (3D), s-wave bound and eontinuum (seattering) solutions of 
Posehl-Teller potentials are obtained by a group theoretie algebraie approaeh 113 ll. 13211 as well as by 
non-algebraie methods 11301.13311. The Sehroedinger equation of relative motion ean be expressed 
in the form 




where 


c/> _ _ __ 

2fi 


7 2 


+ K) 


(9) 


( 10 ) 


By treating the right hand side of Eq. dH) as a souree term, we seek solutions of the homogeneous 
equation 


^r4’±ij') = 0 


( 11 ) 


In the following seetion we present exaet solution of Vo(r). 
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III. SOLUTIONS OF Vo{r) 


To obtain solutions of Eq. (fTTI) . we first convert this equation into the standard equation for 
associated Legendre functions. We then write down the desired solutions as a superposition of 
two linearly independent associated Legendre functions by imposing the boundary conditions for 
scattering states. Next, we rewrite the solutions in terms of hypergeometric functions for asymp¬ 
totic analysis. We then construct Green function of Eq. (fTTI) . In the following subsection we 
briefly describe our method of solution. In passing it is worth noting that as the range ro goes to 
zero, Poschl-Teller potential in ID reduces to a delta well potential Odil . So, it is expected that for 
s-wave interactions three dimensional Poschl-Teller potential will behave as a contact potential. 


A. Scattering solutions of Vq 


Let K = (3/tq, Kr = p, cos 9 = tanh p = z, then we have 


n? 

= -- 


d? cos 6 d 


where 


2/i [dO"^ sin 9 d9 sin^ 0 


q = ik/n = \/—2pE/ k. 


sin^ 9 + Vo 


Lurthermore, substituting Vq given by Eq. (|7]) in Eq. (fTTI) . we have 


(1 - z'^)— - + 1 ) 

dz 


r 


i!{z) = 0 . 


( 12 ) 


(13) 


(14) 


dz"^ dz 1 — z^ 

where the parameter A is given by A(A -I- 1) = 2q!“^ This is the familiar equation of associated 
Legendre functions. Two linearly independent basic solutions of this equation are P^(z) and 

Ql(z). 

By doing asymptotic analysis of the two basic solutions as given in the appendix, one can write 
down the regular scattering solution of the potential Vq in the form 


l/;(z) = Eyx,q 


(15) 


r (i - - - 2) 

with = A~\ where Ag^x is given in the appendix. One can notice that a is bounded by 
l<a<2or0<a<l depending on whether is negative or positive, respectively. In the 
limits tts —)■ ±cxo, a ^ 1 and therefore A assumes integer values of either -2 or 1. 
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B. Green function 


Although Green funetions for one-dimensional Posehl-Teller potentials are studied before OSll . 
to the best of our knowledge the green funetions of Posehl-Teller potentials in three dimension 
(3D) is not eonsidered before. We first note down the Wronskian 113611 




r(i + ^)r(i + ^) 


(16) 


(i-^2)r(i + ^)r(| + ^) 

between P^iz) and QKz). To eonstruet the Green funetion, we prepare another state (p{z) by 
linear superposition between P^{z) and Qx{z) sueh that this state satisfies the outgoing boundary 
eondhion 0(r —)■ oo) I'N-' — gi(fcr+<5o)_ leads to 


r(g) [aA(g) - axi-q)] 


where 


(17) 


«a( 7) = cos[7r(A -f q)/2] sin[7r(g - A)/2] 


(18) 


The Green funetion ean be written in the form 


Ge{p,p') = Ge{z,z') = -TT'l/j{z<)(j){Zy) 


(19) 


That the Green funetion Ge{p, p') is the solution of the Green equation 

■^pGe{p,p) = - ^ - 6{z-z) = p) 


ean be aseertained from the relation 


>'- 1 - 0 + 


' p'—0+ 


dp^pGEip, p') — 1 


( 20 ) 


( 21 ) 


C. Bound state solutions of Vq 

A bound state oeeurs for energy E = —|i7| or A; = i\k\, and so the parameter q defined in Eq. 
(fTSl) beeomes q = —k = —\k\/K. For the solution'0^(r) to behave as abound state, asymptotieally 
^p{r —)■ oo) ~ For negative energy, the asymptotie form of takes the form 

V^(r ^ oo) ~ [Gx{q)e-\'^\^ + Gx{-q)e\^\^] (22) 
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FIG. 1: Absolute square of energy-normalized continuum wave function in unit of is plotted 

against r (in unit of ro) for as = 50ro (solid line), as = lOOro (dashed lines), as = —50ro (dashed-dotted) 
and as = — lOOro (dotted) for A: = 0.01r(^^. 


Now, for r/) to qualify as a bound state, the coefficient of must vanish. Thus, the bound state 
condition is given by 


GxCk) 


2~^T{k) 


t+| r 1+f+f 


= 0 


(23) 


Now, assuming X = -(f + |), and using the identity r(l — x)T{x)7r/ sin[7ra:], we can rewrite 
the above equation in the form 


Gx{k) 



(24) 


This shows that G\{k) =0 when 


— X — k = 2n 


(25) 


with n being an integer, or the argument of the gamma function in the denominator on the right 
hand side of the Eq. (I24l) is a negative integer. As decreases from +oo to 2ro, A decreases from 
-2 to — cxo. This means that the argument of this gamma function can not be negative since k > 0. 
Therefore, the bound state condition is given by Eq. (I25l) which, for low energy, will be satisfied 
forn = 1. On the other hand, when < 0 and as changes from — oo to small negative value A 
decreases from 1 ranging between 0 < A < 1. This means the above equation can not be fulfilled 
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FIG. 2: The phase shift (5o(A:) in unit of vr radian is plotted against dimensionless momentum Zero (scaled 
by a factor of 10^) for three positive and three negative values of scattering length a* as indicated in the 
legend. The lower dotted, dashed and solid curves correspond to = —10^, = —10^ and a* = —10^, 

respectively; while the upper dotted, dashed and solid curves are plotted for positive of same magnitudes, 
respectively. 


for negative seattering length. These findings are eonsistent with the eonditions for bound states 
obtained earlier by group theoretie approaeh [1320. 


IV, RESULTS AND DISCUSSIONS 

In figure 1, we show the square of energy-normalized seattering or eontinuum wave funetion 
\ip{r)\^ di funetion of r for four values of relatively large seattering lengths at a fixed low energy. 
The eontrasting behavior of the wave funetions for positive and negative seattering lengths for 
r > ro is noteworthy. The intereepts of the slope of the wave funetion at r = ro on the x-axis 
are positive and negative for positive and negative seattering lengths, respeetively. The intereepts 
are nearly equal to the respeetive seattering lengths. Thus the wave funetions behave exaetly like 
those of a 3D square well at low energy. At r = ro, the wave funetions for positive and negative 
seattering lengths bend towards and away from x-axis, respeetively. 

Next we analyze phase shifts at low energy. From the asymptotie analysis of the solution for 
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Vo as given in the appendix, we have ^(r —)■ oo) ~ sin[fcr + (5o] where 


So{k) = —tan ^ 

OO 


—ttA / Tik 

cot —-— tanh I — 


1 

Xk 


n=l ^ 


arctan 


k 


-A + n)K 


+ arctan — -— 


arctan ( — 
nn 


(26) 


From the expression for 5o given in Eq. (|2^ one ean notiee that, for fc 7 ^ 0, as A ^ 1 we get 
(5o —^ 7 r /2 — 0+ while as A approaehes to -2, (5o approaches to 7 r /2 + 0+. From the relation 
as = — limfc_j.otan[(5o(/c)]//c, one finds that, for negative the resonance (a^ —)■ — 00 ) will be 
characterized by (5o 7 r /2 — 0+, that is, by A = 1 while for positive the resonance will 

occur when A = —2. The Eq. (|2^ suggests that in order for 6o{k) < 7 r /2 so that < 0, one 
needs to set A = Ai = —1/2 + y/1 + 8a~^/2. Similarly, for > 0 one has to fix A = A 2 = 
—1/2 — \/r+”8a^/2. This analysis shows that we have 0 < (5o < 7 r /2 and 7 r /2 < (5o < for 
Os < 0 and a* > 0, respectively. For |A| « 1 the phase shift 6o{k) (x k in accordance with 
Wigner’s threshold laws. In the limit a ^ 1, A can assume a value of either -2 or 1. In both limits 
as —)■ —cxo and —)■ +cxo, the parameter A is given by A(A + 1) = 2. 

Figure 2 displays the variation of 5o{k) as a function of kro for different values of a*. To 
remain within the validity regime of the effective range expansion, in this plot we have restricted 
the variation of kro at low values below 10“^. From this figure we notice that for large scattering 
length 5o{k) approaches the resonant value 7 r /2 at finite momentum while in the zero momentum 
limit 5o{k) significantly deviates from 7 r /2 even for large scattering length. This implies that 
the momentum- and range-dependence of interactions become particularly important for resonant 
interactions. 

We next discuss how these potentials will be useful to model tunable MFR in ultracold atoms. 
The s-wave phase shift near MFR is 6o{k) = 5bg + <5r where the background phase shift (5bg can 
be approximated as 5bg — —ka^g with Obg being the background scattering length. Here 6r is the 
resonance phase shift. At low energy, the Feshbach resonance width Tj is given by r 7/2 ~ kabgT 0 
where Tq is a parameter related to the width of zero crossing. Now, assuming \kabg\ « 1 one can 
write 


and 


1 

do 


B^-B 


Aa, 


bg 


To 2flbg 


h 


fia'OgT'o 


(27) 


(28) 
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where Bq is the magnetie field at whieh MFR oeeurs (or equivalently, diverges) at zero en¬ 
ergy. Sinee abgFo > 0, the above equation shows that tq will be positive if Obg is positive and 
Obg > ^Jh/(2^VQ) . If ahg < 0 then tq is negative. Sinee the model potentials of Eqs. Q and dS]) 
are derived assuming ro > 0, these potentials will be useful to model those MFR for whieh ro is 
positive. For example, magnetie Feshbaeh Resonanees observed in ^^^Cs near Bq = —11.7, 547 


and 800 G fulfill the conditions for ro > 0 1137. 


3811 . In general, for a relatively broad Feshbaeh res¬ 


onance or an open-channel dominated Feshbaeh resonance [|8l] with positive Obg for which binding 
energy of closed channel bound state has the universal form 77^ ~ ti?/ (2/ia^), the resonance can 
be described effectively by a single-channel potential using our proposed model potentials. The 
treatment of MFR with our proposed finite-range potentials will give new insight into resonant 
phenomena in ultracold atoms by providing effective range dependence of scattering phase shifts 
and near-zero energy bound states. The analytical solutions found in this paper will be useful 
to calculate photoassociation in the presence of an MFR and thereby to develop quantitative un¬ 


derstanding of the Fano effect in photoassociation ll39h . Furthermore, in case of two-component 


Fermi gas of atoms, the modeling of MFR with the finite-range potentials will facilitate to investi¬ 
gate the hither-to-unexplored nontrivial effects of large effective range on the BCS-BFC crossover 
from Bardeen-Cooper-Schrieffer (BCS) state of Cooper-pairs to Bose-Einstein condensate (BEC) 
of dimers. 

V. CONCLUSIONS AND OUTLOOK 


In conclusion we have demonstrated that the effective range dependence of resonant two-body 
interactions can be taken into account within a class of model potentials constructed under effective 
range expansion using the methods of lost and Kohn. We have shown that these potentials reduce 
to the form of Poschl-Teller potential when scattering length is large. We have presented analytical 
scattering and bound state solutions of the potentials for large scattering lengths and established 
the connections of the nature of the solutions with the sign and strength of the A-parameter of the 
Poschl-Teller potential. These finite-range model potentials will permit us to explore the finite- 
range effects of interactions between atoms in low-dimensional traps. In a recent paper lldOll . using 
these potentials we have numerically studied bound-state properties of two atoms in a quasi-two 
dimensional trap and found significant effects of the range on the bound states. For studies of 
collisional properties of cold atoms near resonances, one can employ more accurate multichan- 
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nel scattering-based eomputational methods where one ean use moleeular potentials. However, 
to gain insight into many-body physios near resonanees, one prefers a single-ohannel simplified 
model potential that ean be Fourier transformed so that one ean oonveniently develop many-body 
treatment of a homogeneous system in momentum spaee. Keeping this in mind, we have shown 
that there exists a olass of simplified finite-range potentials that are well suited for the purpose of 
desoribing resonant interaotions more aoourately with both finite-range and energy-dependenoe of 
the soattering proeesses being taken into aooount. 


Appendix A: Asymptotic analysis 


The regular seattering solution is required to fulfill the boundary eonditions ^(r = 0) = 0 and 
ip{r cxo) ~ sm[kr -f (5o(^)] or equivalently, tjj{z = 0) = 0 and ip{z —)■ 1) ~ sin[fi;p -f (5o(/i;)]- 
This means ^jJ{z) will be given by the superposition of the two basie solutions fulfilling these two 
boundary eonditions. Thus we have 


'ipi^z) ~ r 



P^{z) - 




(Al) 


Expressing the two assoeiated Legendre funetions in terms of hypergeometrie funetions, and after 
a lengthy algebra one obtains 




r(i + | 

X z{\ — z^)~^^'^{z^— \)^F 


r(| +j + f) _ 

-|)r(-f-f)sin [vr(| + |)] 



1 

1 

1 

q 1 

1 

1 


3_ 

\2 

2 

2’ 2 

2 

+ 2’ 

2’ ^2 _ 1 


(A2) 


From this expression it is easy to verify that at r = 0 or z = 0, ^ = 0 and therefore is ^ 
regular solution. 

Now, we will show that satisfies the proper asymptotie boundary eondition ^(r —)> oo) ~ 
sin[/i;r -f So{k)]. Sinee z"^ = tanh^p, = — sinh^p, we ean write z{l — z‘^)~^^‘^{z'^ — 1)^ = 
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gi7rA/2 p l^cosh^ p\ and 


1 

1—1 

1 

<< 

1 

1—1 

1 

1 


3 ^2 

V2 

2 

2’ 2 

2 

^ 2’ 

2’z2_i 


F 


X (sinh^p )2 + 2 2F 


r(|)r(g) 


1 

<< 

1 — 1 

1 

<< 

1 

V2 2 

2’ 

2 2 


r(i-| + |)r(i + t + |) 
1 


) 1 ?) ■ . , 2 
Sinn p 


+ 


r(i)r(-g) 

r(i-|-f)r(i + |-f) 


X (Smh p)= = .F(--- + --- + -;l + ,;-^ 


Thus we have 


(A3) 


V^(r ^ oo) ~ A,,, [Gx{q)F'^'^ + Gx{-q)e-^’^^] 


(A4) 


where 


and 


Aq,A = i 


2i+9e*^<ie*^3A/2^i/2p('3/2)r (1 + A + |) 


sin [vr (I + I) r (i + 4 - I) r (-| - I)] 


GA(g) = 


2 -''r(g) 


(A5) 


(A6) 


r(i-| + |)r(i + f + |) 

The quantity within the third braeket on the right hand side of Eq. (A4) ean be expressed in the 
form |GA(q')| sin [kr + (5o] with (5o = f + 0 where 

2-''r(g) 


= arg 


r(i-| + |)r(i + | + f) 


(A7) 


Substituting q = ikj k and making use of the standard phase relationship involving gamma func¬ 
tions we obtain 


0 = argr(i/c/K) — argr(—A + ifc/fi:) — tan 


-1 


—ttA , (irk^ 
cot-tanh — 

2 V 2 . 


(A8) 


Using the formulas [6.1.27], [6.3.16], [6.3.7] of Ref. [l36h . one can obtain Eq.(26). 
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